YOUNG DIAGRAMS AND 7V-SOLITON SOLUTIONS OF THE 

KP EQUATION 
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Abstract. We consider A/-soliton solutions of the KP equation, 



(— Aut + u xxx + 6uu x ) x + 3u y 



0. 



An A/-soliton solution is a solution u(x,y,t) which has the same set of N 
line soliton solutions in both asymptotics y — > oo and y — > -co. The N- 
soliton solutions include all possible resonant interactions among those line 
solitons. We then classify those A/-soliton solutions by defining a pair of N- 



numbers (n + ,n ) with n^ 



(nf, 



,n%), 



£ {1, . . . , 2N}, which labels 



N line solitons in the solution. The classification is related to the Schubert 
decomposition of the Grassmann manifolds Gr(A r , 2N), where the solution of 
the KP equation is defined as a torus orbit. Then the interaction pattern of 
A/-soliton solution can be described by the pair of Young diagrams associated 
with (n+,n _ ). We also show that A/-soliton solutions of the KdV equation 
obtained by the constraint du/dy = cannot have resonant interaction. 



1. Introduction 
In this paper, we consider a family of exact solutions of the KP equation, 

du\ 9 2 "' 



d ( du d 3 u . . 

dx V dt dx 3 dx ) dy 2 



0. 



The solution u(x,y,t) is obtained from the r-function r(a;,y,t) as 

d 2 
u(x,y,t) = 2— logr(x,y,t). 

It is well known that some solutions can be obtained by the Wronskian form (for 
example see 0), 
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where f> n — d n fi/dx n , and {fi(x,y,i)\i ~ 1, . . . , N } is a linearly independent 
set of N solutions of the equations, 

dfj _ d 2 fj dfj _ d 3 f t 
dy dx 2 dt dx 3 



Key words and phrases, the KP equation, N-soliton solution, tau-function, Grassmannian, 
Schubert decomposition, Young diagram. 
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Throughout this paper, we assume fi(x,y,i) to have the form, 

M 

(1.2) fi=^2 °" e<>J ' fOT * = 1 ' • • • ' N > and M > N ' 

with some constants ay which define the N x M matrix Ai N ^ M \ := (ay), and the 
phase functions 6j are given by 

(1.3) e j {x,y,t) = -k j x + k 2 j y-k z j t + e° j for j = l...,M. 

Here kj and 0° are arbitrary constants, and throughout this paper we assume kj 
being ordered as 

k\ < ki < ■ ■ • < ku ■ 

Then choosing particular forms for the matrix -A(jv,m)i one can obtain several exact 
solutions of the KP equation. As the simplest case with N = 1 and M = 2, i.e. 
t = fi = due® 1 + ai2e 02 with anai2 > 0, we have 1-soliton solution, 

(1.4) u = 2^ l °Sr = \{ki~ fc 2 ) 2 sech 2 i(0 1 - 6 2 ) . 

Here we assumed an, a\2 > and absorbed them to the constants #!-, so that we 
have an = a\2 — 1. Note here that if <Xnai2 < 0, then the r-function has zeros and 
the solution blows up at some points in the x-y plane. In the x-y plane, 1-soliton 
solution describes a plane wave u = $>(k x x + k y y — cut) having the wavenumber 
vector k = (k x , k y ) and the frequency u, 

k = (-fci + k 2 , k\ - k\) u = k\-k\. 

Here (k, u>) satisfies the dispersion relation, Auk x + k\ + 3fc^ = 0. We refer to the 
1-soliton solution as a line soliton, which is expressed by a line 0\ = 82 in the x-y 
plane for a fixed t. Then the slope of the line c = dx/dy = —k y /k x represnts 
the velocity of the line soliton in the x-direction with respect to y; that is, c = 
indicates the direction of the positive y-axis. As in the paper [Q, we also call a 
line soliton as [i, j]-soliton, if the soliton is parametrized by the pair {ki, kj). In the 
case of <|1.4ll . the soliton is [l,2]-soliton with the velocity c = k\ + k%. Then the 
ordinary A^-soliton solution consisting of N line solitons, [2j — l,2j']-soliton with 
j = 1, . . . , N, is obtained from the following matrix A^m) with M = 2N 0, 

/l 1 0\ 

1 1 

^4(iV,2A0 = 

\o 1 l) 

Note here that all nonzero entries are normalized to be one by choosing appropriate 
constants (n, and any N x N minor of the matrix Ai N2 n) is non-negative. The 
fact that all TV x TV minors are non-negative (or non-positive) is sufficient for the 
solution to be non singular (see below). 

Also in the paper PP , we found that if all the minors are nonzero and have the 
same sign, the solution u gives an (M — N, A^)-soliton solution consisting of M — TV 
incoming line solitons as y —> — 00 and N outgoing line solitons as y — * 00. The 
set of the r-functions with different size N then provide the solution of the Toda 
lattice. In particular, if M = 27V, we have A^-soliton solution in the sence that the 
solution has the same set of N line solitons in both asymptotics for y — » ±00 (i.e. 
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the sets of incoming and outgoing solitons are the same) . However as mentioned 
in £Q, this iV-soliton solution is different from the ordinary 7V-soliton solution, and 
the interaction of any pair of line solitons is in resonance. In this case, the matrix 
A/n,2N) can be written in the following row reduced echelon form (RREF), 



A(N,2N) 



/l ••• 
1 ••■ 



\0 ••■ !**••■ */ 



where the entries marked by "*" are chosen so that all the N x N minors are 
nonzero and positive (i.e. all the *'s of the last row should be all positive, and 
those in the second row from the bottom are negative and so on). Note here that 
those *'s are all nonzero. Then the line solitons in jV-soliton solution of the Toda 
lattice hierarchy are given by [j, N + j]-soliton for j — 1, . . . , N. 

We then expect that the general case of iV-soliton solutions has a mixed pattern 
consisting of resonant and non-resonant interactions among those line solitons. In 
this paper, we classfy all the possible iV-soliton solutions obtaind by the T-function 
(|l.l(l with fj's in (|1.2(l and M — 2N. The classification is to determine all possi- 
ble patterns of interactions, and it is done by constructing a specific form of the 
coefficient matrix ^4^ :— A( N2 n) m tEH f° r each class of TV-soliton solutions. It 
turns out that a complete classification of iV-soliton solutions can be obtained by 
the complementary pair of ordered sets of iV-numbers, n + = (n^[,...,n N ) and 
n = (nj~, . . . , n N ), which are related to the dominant exponentials in the r- 
function for x — * ±oo and each soliton is given by [n+,n^]-soliton. Notice that 
with the set {1, . . . , 27V}, there are (2iV — 1)!! number of different sets of TV pairs. 
We then claim that each set of TV pairs corresponds to an 7V-soliton solution, and 
provide an explicit way to construct all those iV-soliton solutions. 

It is also well-known that a solution of the KP hierarchy is a GL(oo)-orbit on 
infinite dimensional Grassmannian (Sato Grassmannian) |H]. Then the iV-soliton 
solution can be identified as the 27V-dimensional torus orbit on the Grassmannian 
Gr(iV, 2iV), and the coefficient matrix An represnts a point on Gr(JV, 2iV). Since the 
Grassmannian has the Schubert cell decomposition, one can first classify the matrix 
An using the decomposition, i.e. identify the cell An belongs to. This classification 
corresponds to the set n + , and the n~ gives a further decomposition of the cells 
to classify the interaction patterns of iV-soliton solutions. We then define a pair of 
Young diagrams {Y + , Y~) associated with the pair of the number sets (n + , n~). In 
particular, the number of resonant vertices in the iV-soliton solution parametrized 
by [Y + , Y~) is given by 

2 I I I I ' 

where |K | denote the size (degree) of the diagrams, i.e. the number of boxes in 
the diagrams ( Theorem 14. 3|) . 

The paper is organized as follows: In Section 2, we give a general basic structure 
of the T-function, and present all types of 2-soliton solutions, which provide the 
building blocks for the general case of 7V-soliton solutions. Here we also introduce 
the pair of numbers (n + , n~), and the Young diagrams (Y + , Y~). In Section 3, we 
briefly introduce the Grassmannian Gr(7V, M) and the Schubert decomposition of 
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Gi(N, M). We also identify iV-soliton solution as an 27V-dimensional torus orbit on 
Gr(7V, 2N), and briefly mention that Gr(N, 2N) contains all possible (m, n)-soliton 
solutions for 1 < m < N and 1 < n < N, which are distinguished by the Schubert 
decomposition. Here (m, n)-soliton is the solution consisting of m incoming solitons 
for y — > — oo and n outgoing solitons for y — > oo (see £Q). In Section 4, we describe 
the structure of the coefficient matrix An for each ./V-soliton solution by prescribing 
an explicit construction of the matrix An- We here define N-soliton condition on 
the matrix An which determines local structures based on the types of 2-soliton 
interactions in the ./V-soliton solution. Finally, in Section [5J we discuss iV-soliton 
solutions of the KdV equation, and show that the KdV iV-soliton solution cannot 
have resonant interactions. 



2. Basic structure of the t-function and 2-soliton solutions 

Let us start with the follwoing Lemma which shows the basic structure of the 
T-function: 



Lemma 2.1. The T-function \l.l\l with fj 's given in jj.jjj) can be expanded as 
sum of exponential functions, 

t= ^2 ^(ii, - - - , ijv) JJ (fctj-feijexp 

l<ii<— <i N <M \<j<l<N 

where £(zi, . . . , «at) is the N x N minor given by the ij-th columns with j 
in the matrix A^ n ,m) — i a ij) °f HIW > 




.N 



€(*!,-- • ,«jv) := 



°i,»i 



0,1,: 



OjV,ii ' ' ' 0,N,i N 

Proof Apply the Binct-Cauchy theorem for the expression (see for example : 3: ), 
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where E) ' = (-kj^Ej with E, 

From this Lemma, it is clear that if all the minors £(ii, . . . , «at) are non-negative 
(or non-positive), then the r- function is sign definite and the corresponding solution 
u is nonsingular (recall also the order k\ < ■ • • < fejw)- This is the first restriction on 



the matrix A 



(N,M) 



for 7V-soliton solution, and we later impose that all the N x N 



minors are non-negative after making Ar N .M) to be in the row reduced echelon 
form. 

We now consider the case of ./V-soliton solutions, that is, we take M = 2N. 
With the ordering of the numbers fc/s, i.e. fci < fc 2 < • ■ • < k 2 N, one can find the 
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asymptotic behavious of the r-function for x — ► ±00: For this purpose, we first 
define the function w given by 

d 
w(x,y,t) = -— logr(x,y,t). 

ox 

Suppose that w has the asymptotic form with some numbers {n 1 , . 
{!,..., 27V}, 



'A' 



} c 



(2.1) 



N 



AT 

.7=1 



as 



as 



.x — * —00 . 



It is clear that nf ^ n+ and n i ^ n- if i ^ j. Then with the sets of iV numbers 
s in the asymptotics i|2.1[l , we define the following two N- vectors with the entries 



,±> 



from the set {1, 2, . . . , 2N}: 

Definition 2.1. A pair (n + ,n 
n+ = (n+,n+,.. ■■' 



K , n 2 , . 



AO 



l NJ 



of iV- vectors are defined by 
with 1 = n~l < n\ < ■ ■ ■ < n^ < 2N, 



with n • > nj . 



The ordering in the set n + is just a convenience, and once we made the ordering in 
n + , we do not assume a further ordering for other set n~. We assume that those 
sets are complementary in the set {1,2,..., 2iV}. As will be explained below, this 
assumption is necessary to have an Af-soliton solution, that is, the solution has the 
same sets of N line solitons for both asymptotics y — > ±00. 

With those number sets n + and n _ , one can define the pair of Young diagrams 

(Y+,Y-y. 

Definition 2.2. We define a pair (Y + ,Y~) associated to the pair (n + ,n~): 



The Y + is a Young diagram represented by 



£,. 



,£+) with ef = nf-j 



where each I , represents the number of boxes in a row counting from the 



bottom (note £+ < if +1 )- 
• The Y~ is a Young diagram associated to (t[ 
as 



£ N ), where t, is defined 



■ ,^-jy 



< n k n k < n- for k > j > 



Then I ■ represents the number of boxes in a row arranged in the increasing 
order from the bottom. 

Notice that the Young diagram Y + is uniquely determined from n + , but the 
diagram Y~ is not unique (note that (£j~, . . . ,£jf) is unique of course). However 
the number of boxes in Y~ gives the number of intersection vertices of iV-soliton 
solution having a particular type of interaction ( Theorem 14.3(1 . 
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As a main goal of this paper, we will show in Section^that for each pair (n + , n~) 
one can construct an AT-soliton solution consisting of N line solitons, each of which 
is given by [n^~,n^]-soliton for j = 1, . . . ,N. The total number of different N- 
soliton solutions is given by (2N — 1)!! which is the number of different sets of 
pairs { (n + , n~) }. Namely, we construct an ./V-soliton solution which is labeled by 
(n + ,n~) as follows: 

1) Consider a set of 2N numbers {1, . . . ,27V}, which are associated to the 
parameters k\ < ■ ■ ■ < k,2N- 



Choose N different pairs from {!,..., 27V}, and label each pair as 



n h n j 



so that the numbers n- satisfy 
i) nf < nj for j = 1,...,7V, 
ii) 1 = nf < nf < • • • < riff. 
3) Define the pair (n + , n~) with n^ = (nf, . . . , nf), and construct the corre- 
sponding matrix An — A/n,2N)- 

Note here that the labeling is unique with a given choice of pairs, and the corre- 
sponding iV-soliton solution consists of [nf, nJ]-solitons for j = 1, . . . , N. The item 
3) will be given in Section 0] 

In order to classify those iV-soliton solutions, we first note: 

Lemma 2.2. For each given ordered set n + , the number of choices o/n~ is given 

by 

N 

m n+ =l[(2j-n+). 

3 = 1 

Proof. Since n~^ > n N , there are (2iV — n N ) numbers of possible choices of n N . 
Having made a choice of n N , one has 2(7V— 1) — n N _ l many choices for n N _ 1 . Now 
repeating this, the result is obvious. □ 

This Lemma provides the number of different 7V-soliton solutions having the 
same n + , that is, the functions w — —(d/dx)\ogr for those solutions have the 
same asymptotic values for x — > ±oo. As a corollary of this Lemma, we also have 
the following identity, 

(2JV-l)!! = ^m n+ . 

We also have the following Lemma for the ordered set n + : 

Lemma 2.3. Each number nj in the n + is limited as 

n+<2j-l for j = l,...,N. 

Proof. Since nf < nj , there are N — j numbers of nf and N — j + 1 of n^ , 
which are larger than nf , that is, we have 

(N - j) + (N-j + l)<2N- nf . 

This implies the Lemma. □ 

In terms of the Young diagram Y + , we have 
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Corollary 2.1. The maximum diagram for Y + , denoted as 7 max , is the Young 

diagram associated to n + = (1,3, .. . ,27V — 1), that is, Y max is the up side down 

staircase with the size, 

N(N-l) 

1 max | — ~ 

The maximum Young diagram for Y~ is also given by Y max , which now corresponds 
to the case n = (27V, 27V — 1, . . . , 7V+ 1). Then the Young diagrams Y ± associated 
to N -soliton solutions are given by a subdiagram ofY max . 

Note that |Ym ax | is a l so the number of intersection vertices of TV lines in a plane 
in the general position. We will then determine whether each point is of resonant 
or non- resonant. 

In Figure ^ we illustrate the asymptotic stage of an 7V-soliton solution with TV 
solitons marked by [rif, n~] for j = 1, . . . , N. We note a duality of the values of w in 

the sense that if w takes Y] - —1 fej. , denoted as (ii, . . . , «at) in the Figure, then w also 
takes (ji, . . . , jn) which is the complementary set of (ii,.. ■ , iw) in {1, ... , 27V}. In 
terms of the minors, this implies that if £(ii, . . . ,«tv) ^ then £(ji, . . . , jn) i= 0. 
The duality can be considered as a symmetry of the KP equation: In particular, if 
9® = for all j, then TV-soliton solution has the symmetry (x, y, t) ^ {—x, —y, —t). 
In this case, at t = all the solitons intersect at the origin in the x-y plane (see 
Figure |2J. Then each soliton can be moved by shifting 6^. With this symmetry, we 
define the duality of the r-function: 

Definition 2.3. We say that the r-function in Lemma \l. II with M — 27V satisfies 
the duality, if for any complementary sets {i±, . . . , ij\f} and {ji, . . . , j/v}, the minors 
satisfy 

£(«i,-..,*jv) =0 if and only if £(ji, . . . ,j N ) = . 

Before we discuss the general case with arbitrary TV, we list up all possible 
2-soliton solutions, and show that the r-functions for those cases all satisfy the 
duality. As will be shown below, those solutions give the building blocks of 7V- 
soliton solutions. There are three cases of the solutions (see also p^), which are 
labeled by the pair (n + , n~) with n ± = (nf , n^ )■ Two line solitons are labeled by 
\n\,n{\ and [n+,n^]: 

i) n + = (1,3) and n~ = (2,4): This corresponds to the ordinary 2-soliton 
solution, and the matrix Ai which we denote as A 2 , takes the form (after 
choosing 6°: properly), 

jo) _ (I 1 

2 " [p 1 1 

which gives only 4 nonzero minors, £(1,3), £(1,4), £(2,3) and £(2,4). This 
shows the duality, i.e. £(1,2) = £(3,4) = 0. Those nonzero minors then 
label the four asymptotic regions separated by two line solitons: Denoting 
(i, j) := ki + kj, w takes (1, 3) as x — > oo which is dual to (2, 4) as x — v — oo. 
Also w takes (1, 4) as y — ► oo, which is dual to (2, 3) as y — > — oo. This also 
implies that there is no other transition of the w value, hence the interaction 
is not a resonant type. 
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Figure 1 . An 7V-soliton solution. The label [nf , n, ] indicates 
the [n+,n~]-soliton. Here the soliton velocities Cj := k n + + k - 

3 3 

are assumed to be ordered as Cj+i > c,j. Each asymptotic region 
is marked by the value of the function w there, i.e. (i%, . . . ,in) 
implies w = J2j=i *V 



The Young diagram Y + in this case is the maximum one, i.e. Y + = 
□. The diagram Y~ has no box, Y~ = 0. The interaction of O-type is 
thus labeled by the Youg diagram Y + . It is also useful to note that the 
intervals defined by the labels [1,2] and [3,4] of solitons have no overlaps. 
Non-overlapping property can be stated as £(1,2) = and its dual one 
£(3,4) = 0. 
n + = (1,2) and n = (3,4): This case corresponds to the 2-soliton solu- 



tion of the Toda lattice , and the matrix Ai , denoted as A 2 
structure, 

'1 

,0 1 

where "+ 



(T) 



has the 



4 T) = 



example of A 2 is 



shows the signs of the entries (also nonzero). An explicit 



4 T) = 



and the corresponding 2-soliton solution is illustrated as T-type in Figure 
|5J Following the arguments in [J, one can easily show that two solitons 
are given by [1,3]- and [2, 4]-solitons, and the interaction is in resonance. 
The main point in this example is that all the 2x2 minors of A 2 ' are 
nonzero, and this results resonance interaction. Non-vanishing condition 
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O-type 





P-type 



Figure 2. 2-soliton solutions. Two line solitons in those figures 
are [1, 2]- and [3, 4]-solitons for O-type, [1, 3]- and [2, 4]-solitons for 
T-type, and [l,4]-and [2, 3] -solitons for P-type. The parameters 
are chosen as (ki, fe2, fe, ki) = (—2,0,1,2) and 6° = 0, Vj for all 
types. 



of all the minors is necessary to make the resonsnt interaction, and those 
minors describe the intermediate solitons which forms Y"-shape resonant 
interaction (see pQ for more detail). 

Again we notice the duality in the values of w in the asymptotics: w 
takes (3,4) for x — + — oo which is dual to (1,2) for x — > oo. Also w takes 
(1,4) for y — ■> oo, which is dual to (2,3) for y — > — oo. One more duality 
appears in the inside of the resonant hole, that is, w takes (1, 3) for t > 0, 
which is dual to (2,4) for t < (or vice versa). Thus all six nonzero 
minors contribute to make the resonant interaction. The pair (Y + ,Y~) of 
the Young diagrams is then given by (0, 0). Alos note that the intervals 
defined by the labels of solitons, [1,3] and [2,4] have a partial overlap, 
iii) n + = (1,2) and n~ = (4,3): This case has been noted in jH], and the 
corresponding matrix A 2 , denoted as A^ ("P" stands for Physical, see 
Remark 12.411 . is given by (after choosing 9® properly), 



A 



(P) 



10 0-1 
110 



which gives again only four nonzero minors with £(1, 2), £(1, 3), £(2, 4) and 
£(3, 4). This implies also no resonance, and the solitons are given by [1, 4]- 
and [2, 3]-solitons. This solution is similar to the O-type, but the pair of the 
Young diagrams is different, and given by (0,D). This interaction of this 
type is labeled by the Young diagram Y~. Also the intervals defined by 
the labels of solitons, [1,4] and [2,3], have complete overlap. This implies 
£(2, 3) = and its dual f(l, 4) = 0. 

One should note that we have the following order in the soliton velocities Cy := 
ki+kj among those solitons for given constants kj with the order ki < k 2 < k 3 < k 4 , 

Cl2 < C13 < C14, C 23 < C 2 4 < C34 . 

Here note that Cu and C23 cannot have a definite order (even cu = C23 is possible) . 
This ordering indicates that any 2-soliton solution can be classified as one of those 
three solitons. One should note here that the type for a given pair of solitons are 



completely determined by their labels 



"] for j = 1,2. Namely the inetraction 
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pattern is completely determined their asymptotics in y — » oo, where two solitons 
are labeled. 

Although the interaction patterns of the cases of O-type in ii) and P-type in iii) 
are the same, i.e. nonresonant, they describe orbits on different cells in terms of the 
Schubert decomposition of the Grassmannian, in which each cell is parametrized 
by a Young diagram. We will explain more details on this description in Section 

Remark 2.4. In [7], the KP equation was introduced to describe a transversal sta- 
bility of the KdV soliton propagating along the ir-axis. Then the order of the 
wavenumber in the y-direction is assumed to be much smaller than that in the 
x-direction. This implies that the wave described by the KP equation should be 
almost parallel to the y-axis for a better approximation within the physical setting. 
In this sense, the P-type solitons are more physical than other two types. The 
KdV two solitons are also obtained from the P-type solitons with k\ < hz < and 
fo = —k 2 , k^ — —k\ (see Section EJ. 

3. Grassmannian Gr(A, 2A) and A-soliton solutions 

Here we briefly summarize the basics of the Grassmann manifold Gr(A, M) in 
order to explain that the minors £(?i, . . . , in) hi Lemma |2. II provide a coordinate 
system for Gr(A, M), the Plucker coordinates. Namely a solution given by the r- 
function i|l.l|) can be marked by a point on Gr(iV, M). The purpose of this section 
is to identify an A-soliton solution as a 2A-dimensional torus orbit of a point on 
Gr(A, 2A) marked by the minors £(?i, . . . , In)- Then using the Schubert decompo- 
sition of the Grassmannian Gr(A, 2AT), we classify the orbits which represent the 
A-soliton solutions. 

3.1. Grassmannian Gr(N,M). A real Grassmannian Gr(A, M) is the set of TV- 
dimensional subspaces of R M . A point £ of the Grassmannian is expressed by the 
A-frame of vectors, 

M 

i = [&,&, • • • , £n], with & = J2 a ljSj G K M , 

where {e, | i = 1, 2, • • • , M} is the standard basis of R M , and (a^) = Ar NiM \ is the 
N x M matrix given in (|1.2|) . Then the Grassmannian Gr(A, M ) can be embedded 
to the projectivization of the exterior space /\ M M , which is called the Plucker 
embedding, 

Gr(A,M) <-» P(A W IR M ) 

Here the element on P(/\ K M ) is expressed as 

CiA---A^tv= ^2 €(h,--- ,ijv)e n A ••• Ae lN , 

l<ii<---<i N <M 

where the coefficients £(ii, • • • , in) is an N x A minors given in Lemma l2.ll which 
are called the Plucker coordinates. 
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It is also well known that the Grassmannian can have the cellular decomposition, 
called the Schubert decomposition ^Q, 



(3.1) Gt(N,M) 

where the cells are defined by 



u 



W(ii,--- An) 



W( 



n, 



Kil<---<ijv<M 



/0 



Jn) = { € = 



°\ 









• 


■ 


1 





• 


• 


* 





• 


• 





1 


• 


• 


* 


* 




• 








• 


• 1 



A' 



.M 



t,M 



V 



{ the matrix £ = AT N M s in the echelon form 



whose pivot ones are at (i\, ■ ■ ■ , ijy) positions } 
Namely an element £ = [£i, • • ■ , £,n] £ W(ii, ■ • ■ , iff) is described by 



f eW(»i 



,i N ) •» 



(i) H(i u --- ,i N )^0. 

(ii) £(ji, • • ' ,3n) = if j„ < i„ for some n. 



Each cell W{i\, ■ • ■ , «at) is called the Schubert cell, and it is convenient to label the 
cell by the Young diagram Y = (i\, ■ ■ ■ ,«&) where the number of boxes are given 



by 



] for j = 1, 



, N (counted from the bottom) , which also expresses 



^N 



a partition (£n,£n-i, ••■ ,£i) of the number \Y\ := J2i=i ^i> tn e size of Y. The 
codimension of the cell W(ii, ■ • ■ , ife) is then given by the size of Y, i.e. 

codimW(ii,--- ,i k ) = \Y\, 

and the dimension is given by the number of free variables in the echelon form. 
Note that the top cell of Gr(7V, M) is labeled by Y = (1, • • ■ , N), i.e. \Y\ = 0, and 

dim W(l, ■■■ ,7V) = dim Gr(/V, M) = N x (M - N). 

The main point here is that each matrix Ai NM \ in (|1.2|) can be identified as 
a point on a Schubert cell of the Grassmannian. Now we identify an iV-soliton 
solution as a torus orbit of a point on the Grassmannian Gr(iV, 27V), and show that 
orbits in the different cells give different types of 7V-soliton solutions. 

3.2. 7V-soliton solution as a 27V-dimensional torus orbit. First we note that 
the r-function i|l.lfl in the Wronskian form can be expressed as 

(3.2) t(x, y,t) = ( v N , KD{x 7 y, t)A* ■ v N ) , 

where jj" is the highest weight vector in a fundamental representation of GL(27V) 
which is given by an element of /\ M. 2N , 

v = e\ A 62 A • • ■ A ejv , 



12 YUJI KODAMA 

and K, A (A* represents the transpose of A) are 2N x 2N constant matrices denned 
by 

K — (kij)i<i,j<2N , with kij := (-fcj) 4 " 1 

A = (aij)i<i t j<2N ■ 
The D(x,y,t) is 2N x 2N diagonal matrix associated to an (R*) 27V -torus action, 

D(x,y,t) = diag(e e \...,e e2JV ) with 9j = -kjX + k)y - k)t + tf . 

Note here that A 1 ■ v N defines the matrix ^4^ = A( N2 n) m JT75J, that is, 

A f ■ ei A • • • A e N = A*ei A A'e 2 A . . . A'e^v , 

which can be idetified as ^4^,. The inner product (•, ■} on /\ M. 2N is defined by 

(vi A • • • A Vj, Wi A • • • A lu,-) = det [((u m , tw n ))i<m,n<i] , 

where (v m ,w n ) is the standard inner product of v m , w n G M. 2N . 

The expression (|3.2() implies that the solution u(x,y,t) given by the r-function 
(|l.l|l is a 2iV-dimensional torus orbit, (M*) 2Ar -orbit, of a point on Gr(iV, 2N) marked 
by An- Note also that the Schubert cell W(i\,... ,ijy) is invariant under this 
torus action. The orbits, iV-soliton solutions, can be first classified in terms of the 
Schubert decomposition with the index sets (ii, .. . ,i/\r). 

Example 3.1. Gr(2,4): The Schubert decomposition of Gr(2,4) is given by 

Gr(2,4)= [J W(i,j). 

l<ij<4 

There are six cells W(i,j) with dimW(i, j) = 7— (i +j). 

i) W{1, 2): This is a top cell of maximum dimension four, and a point on the 
cell is given by 

(\ * *\ 
2 ~ \0 1 * *) ' 

which includes 2-soliton solutions of T-type, A 2 , and P-type, A 2 ' . 
ii) W(l, 3): A point on this cell is marked by 

(\ * *\ 
2 ~ ^0 1 *) ' 

The 2-soliton solution of O-type with A 2 belongs to this cell. Thus O- 
and P-types belong different cells, even though they show nonresonant in- 
teractions, 
iii) W(l,4): The A 2 matrix has the form, 

(\ * * o\ 
2 _ ^0 l) ' 

In this case, nonzero minors are given by £(i, 4) for i = 1,2,3. This im- 
plies that the exponential term E4 = e 4 becomes a common factor in the 
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r-function. Then the solution gives (2, l)-soliton solution, that is, two in- 
coming solitons and one outgoing soliton in the asymptotics y — > ±00 (see 

00) 

iv) W(2, 3): The A 2 matrix is given by 



An 



Then the exponential term E\ = e 01 is missing in the r-function, and the 
solution describes a (1, 2)-soliton solution. 
v) W(2,4): The A 2 matrix is 

A 2 =(° l * °V 

^0 1) 

The solution is a 1-soliton solution of [2,3]-soliton. 
vi) W(3,4): The A 2 matrix is 

A 2 =(° ° l °). 

\o 1) 

The r-function is just a product of E 3 = e 03 and E4 = e 9i , and the solution 
is trivial, u = 0. 

Remark 3.2. Since W(l,3) is a boundary of the top cell W(l, 2), a point on PF(1,3) 
can be obtained a limit of a sequence of points in W(l, 2). For example, consider 
the following matrix (representing one parameter family of points in W(l, 2)), 

4>(e) = 

where the parameter e is a positive number. If e < 1, then the matrix gives 2- 
soliton solution of T-type (note that all the minors are positive) . Then in the limit 
e — > 0, the matrix can be interpreted as the O-type matrix A 2 in W^(l,3), i.e. 
(recall here that e~V = r is the Pliicker coordinate, i.e. they both give the same 
solution). It is also obvious that A 2 is obtaind by a limit of T-type matrix. 
Thus T-type is generic, while the O-type and P-type are non-generic. Also in the 
Gallipoli conference in June-July 2004, Pashaev showed that 2-soliton solution of 
T-type is a degenerate case of 4-soliton solution of nonresonant type (obtained by 
the Hirota method) [Sj. 

Since ./V-soliton solution is an orbit on a particular Schubert cell W(i\, . . . , in) 
which is a boundary of the top cell W(l, . . . , TV), it seems to be true that all the 
iV-soliton solutions are obtained by the limits of the r-function of T-type which 
lives on the top cell. We will discuss more details in a future communication. 
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4. Construction of the An matrix 

In this section, we give an explicit construction of the matrix An for 7V-soliton 
solution parametrized by the pair (n + , n _ ) : The key idea for the construction is 
based on the classification of 2-soliton solutions, that is, a local structure of the 
matrix contains one of those types, 0-, T- or P-types. With a given pair (n + , n _ ), 
we construct an iV-soliton solutions consisting of N line solitons with the labels 



n t> n S 



] for j = 1, . . . , N. In order to do this, let us first define: 



Definition 4.1. The length of [n^,n- ]-soliton, denoted by Lj, is defined by 

Lj :— min { nj — nj , 2N — nj + nj j . 
Note that 1 < Lj < N for all j = 1, • • • , N. 

In the cases of 2-soliton solutions, we have 

i) for O-type, L\ = L% = 1, 
ii) for T-type, L\ = L<x = 2, 
iii) for P-type, L\ = L<i = 1. 

Then we note that the length of the line soliton gives the condition on the minors. 
In those cases, £( n t > n j) = if and only if Lj < 2. 

For iV-soliton solutions with [n^,n^]-solitons, we impose that the r-functions 
satisfy the following condition which plays a crucial role for the construction of Am'- 

Definition 4.2. We say that a r-function satisfies the N-soliton condition, if it 
satisfies the duality (Definition 12.3(1 . and has missing exponential terms (i.e. van- 
ishing minors) if and only if there exists [n^nj] -soliton whose length is less than 
N. If Lj < N for some j, we require that 

• for Lj = nj — Hj, then the TV x N minors satisfy 





£(..., [nj.nj + l,. 


■■,Uj] 


,...) = o, 


• for Lj 


= 2N - nj +n+, then 








£([l,...,n+], 


, [nj,.. 


• , 2N] ) = 



Here the Lj + 1 entries inside the brackets [• • ■ ] are consecutive numbers, and the 
other entries in the ". . . " outside the [• ■ • ] can take any numbers. 

Recall that in the case of 7V-soliton solutions for the Toda lattice, every soliton 
has the length N, and there is no vanishing minor pp. The 7V-soliton condition 
is a natural extension of that in the case of 2-soliton solutions, and it provides a 
local information of 2-soliton interactions in 7V-soliton solution which are classified 
into three types, 0-, T- and P-types. Note in particular, with this condition, one 
identifies all the vanishing minors in the r-function. 

Let us now construct the matrix ^4^ The construction has several steps, and it 
might be better to explain the steps using an explict example. 

Let us consider the case N — 4 with n + = (1, 2, 4, 5) and n~ = (3, 7, 6, 8): 
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Step I: Since the minors £(1,2,4,5) ^ and £(3,6,7,! 
that An matrix has the structure, 

/* 



t^ 0, one can easily see 



(4.1) 






* * * 





0* 


v° 






0\ 



*/ 



where "*" marks a possible nonzero entry, and in particular, the first element in 
each row should be nonzero (this is the pivot place). Note here that £(1, 2, 3, j) = 
for any j, which is also considered as the 4-soliton condition with the length L\ = 2 
of [l,3]-soliton. Notice that £(«i, i 2 , H^^) — for any 4 < i 1; . . . , i 4 < 8 are dual 
to £(1,2,3, j) — for some j. 

Step II: In this step, we put the form l|4.1[l into the RREF. First we normalize 
the first nonzero entries to be 1 in every rows, and eliminate the (2,4)-entry by the 
third row and the (3,5)-entry by the forth row. Note that this process does not 
affect the 0's of the matrix in (|4.1|) . Then we have 



A 



v° 



0\ 



V 



Further eliminating "*" in the 2nd columns by the 2nd row, we have 



(4.2) 



A 





# 
# 
* 



# 
# 
* 



# 
* 



where the two rows in the 
£(5,6,7,8) = £(4,6,7,8) = 



2x3 submatrix with # entries are parallel (recall 
0). Also note that the signs of nonzero entries are 
uniquely determined by the requirement of all the minors to be nonnegative, for 
example, the last three entries in the first row are nonpositive, and the signs of the 
last three entries of the following rows change alternatively. Up to this point, we 
did not use the order in n . In another words, the form of the matrix in 14.2fl is 
common for any n - with a given n + . 

Step III: Now we use the orders in n - to determine an explicit entries in the 
matrix for an 4-soliton solution consisting of four line solitons labeled by \nj , nj], 
i.e. [1,3], [2,7], [4,6], [5,8]. In this step, we identify all other zero elements (i.e. 
zero minors) of the matrix An using the iV-soliton condition fDefinition l4.2|l . The 
follwoing graph shows the overlapping relations among those solitons and the length 
of each soliton: 



(4.3) 



o — — 



Li = 2 
L 2 =3 
L 3 = 2 
L 4 = 3 



Now we use the 4-soliton condition. First we notice that from L\ = 2 < 4 we have 
£(1,2,3, j) = for any j (see Step I). Also from L$ — 2, we have £(4, 5, 6,j) = for 
any j. In particular, £(1, 4, 5, 6) = implies that (2, 6)-entry must be zero, so that 
(l,6)-entry also zero. Also £(3,4,5,6) = gives £(1,2,7,8) = which is also the 
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FIGURE 3. 4-soliton solutions labeled by n+ = (1,2,4,6) and 
n = (3,7,6,8). The left figure shows the 4-soliton solution 
at t = and the right one at t — 10. The parameters 
are given by 9® = for all j = 1, . . . , 8, and (fci, . . . , k&) = 
(-3/2,-1,-1/2,0,1/2,1,3/2,2). The solution at t = -10 is ob- 
taind by the reflection (—X, —y) of the right graph. 



4-soliton condition with L 2 — 3. The £(1, 2, 7, 8) = implies the 2x2 determinant 
of the right bottom corner must be zero. We then get the following structure, 

(l * # jf\ 

1 * # # 

1 * x x 

\0 1 * x x/ 

where the 2x2 minor with the x entries is zero. Also from the partial overlaps (i.e. 
T-type interaction) in [1,3] with [2,7], in [2,7] with [5,8] and in [4,6] with [5,8], 
we need nonzero entries for those marked by *, =ff and x. Then one can easily find 
an explicit example of the matrix A^ whose 4x4 minors are all nonnegative, 



(I 





-2 











-3 


-9\ 





1 


1 











1 


3 











1 





-1 


-2 


-2 


1° 











1 


1 


1 


1/ 



In Figure |21 we illustrate a 4-soliton solution generated with this matrix. Notice 
that there are three resonant vertices corresponding to the interactions between 
[1, 3] and [2, 7], [2, 7] and [5, 8], and [4, 6] and [5, 8] (see the graph[OJ|. Other three 
vertices are two of O-types and one of P-type. 

One can show the following Proposition for the general case: 



Proposition 4.1. // there is a pair [n^,rij ] with Lj = 



< N, then the 



N x (Lj + 1) submatrix (a mn ) with 1 < m < N and nj < n < n- of the matrix 
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An has the maximum rank Lj. (Similar thing is also true for the case Lj 
2N -n~ +nf.) 



Proof. Since Lj < N, we have £([r 



(see the 7V-soliton 



condition). Here the entries marked with * are arbitrary columns of An whose 



maximum rank is TV. If we choose N 



1 independent columns for those 



entries, then the £ = implies that the submatrix can have at most the rank Lj. 
Then the following is imediate as a Corollary of this Proposition: 



D 



Corollary 4.1. If a pair [n'J.nJ] has the minimum length, i.e. nj 
the N x 2 submatrix with the columns (a mn ) with 1 < m < N and 
the matrix should be in the form, 



1 



n = rif, 



then 
i7 of 



"o 


o" 








1 


1 











0_ 



where l's are in the jth row. 



Proof. From Proposition l4.il the rank of the submatrix is one. This implies that 
the rows of the submatrix are all proportional. Since the entry (j,n^) is the pivot 
1, all the other entries of this column should be zero in the RREF. This gives the 
result. □ 

This Corollary shows that the line soliton marked by \nj ,n~] of length one 
intersects with all other solitons without resonance, that is, the interactions are 
either O-type or P-type. Notice that in this case [nj ,nj] cannot have partial 
overlap with other pairs, and this can be seen in the matrix An- Proposition 14. II 
provides the interaction structure of iV-soliton solution for O-type and P-type. 

We now show that the Young diagrams Y + and Y~ provide some information 
on the interaction pattern of 7V-soliton solution. If we ignore the phase shifts and 
resonances among the solitons, 7V-soliton solution has — - 2 ~ ' vertices in a generic 
situation. Then one can find the number of vertices having a particular type of 
interactions. 

Theorem 4.3. The N -soliton solution labeled by the pair of Young diagrams (Y + , 
Y~) has 

• \Y + \ number of vertices of O-type, and 

• \Y~\ number of vertices of P-type. 

The number of resonant vertices, T-type, is then given by 

N(N - 1) 



Y^ 



\Y- 



Proof The Young diagram Y + can be obtained directly from the matrix An'- 
For example, consider n + = (1,2,4,5,7) and n~ = (3,9,6, 10,8). Then the graph 



18 YUJI KODAMA 

showing the overlaps among [ru, nj] is given by 
o — • 

o — — — — — — • 

□ o — • 

□ o — — — — • 

□ □ o • 

where the boxes under the non-pivots (»'s) give Y + , which is given by (1,1,2) 
representing the number of boxes in the row from the bottom (Y + shape is just 
up-side-down picture of one given in the graph). It is then obvious that each □ 
provides O-type interaction, i.e. two solitons having no overlap. For example, [1,3] 
has no overlap with [4,6], [5, 10] and [7,8]. This proves that \Y + \ gives the number 
of vertices of O-type. 

The diagram Y~ gives the number of reverce in the sequence (n7/ , . . . , njf). Again 
consider the example above. The number n^ = 9 has two reverces, 6 and 8. Then 
from the graph above, this implies that [2,9] has complete overlap with [4,6] and 
[7, 8], i.e. [2, 9]-soliton has P-type interaction with those solitons. Now it is obvious 
that \Y~\ gives the number of vertices of P-type. □ 

Before ending this section, we give a complete list of the 3-soliton solutions with 
explicit examples of A3 matrices: 

Case a: n + = (1, 2, 3): Then we have m n + = 1-2-3 = 6 different choice of n~. 
that is, six different 3-soliton solutions. Since the Young diagram Y + has no box, 
Y + = 0, there is no O-type vertex in those solutions. 

al. n = (4,5,6): This is the case of Toda lattice pQ. The pair of Young 
diagrams is (0,0), that is, all the vertices are T-type. Example of the 
matrix A3 is given by the following matrix. We also show the overlapping 
graph. 








1 


3 


6 \ 





— 


— 


• 







1 


-1 


-2 


3 , 







— 


— 


• 





1 


1 


1 


1 









— 


— 



a2. n = (4,6,5): This has two T-type vertices and one P-type vertex, i.e. 
Y~ = □. In particular, the P-type interaction between [2, 6]- and [3, 5]- 
solitons implies that the (3,6) entry should be zero (because of £(1, 2, 6) = 
0). Then we have an example 











1 


2 


4 \ 





- 


- 


• 







1 





-1 


-2 


-3 , 







— 


— 


— 








1 


1 


1 












— 


• 



Note here that the duality of £(1, 2, 6) = and £(3, 4, 5) = 0. Both vanishing 
minors are due to the 3-soliton condition for L% = L3 = 2. 
a3. n = (5,4,6): This has again two T-types vertices and one P-type vertex 
(notice that the Young diagram Y~ is the same as the case 2a). An example 
of A3 is 



1 








1 


2 \ 





— 


— 


— 


• 





1 


-1 


-2 


- 2 ■ 







— 


• 







1 


1 


1 


1 









— 


— 
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Again notice the 3-soliton conditions for L\ = L2 = 2 which give £(1, 5, 6) = 
and its dual £(2, 3, 4) = 0. The £(2, 3, 4) then gives (1, 4)-entry to be zero. 
a4. n = (5,6,4): There are two P-type and one T-type vertices. We have 

'1000 1 1 \ o - - - • 

10 0-2-1, o - - - • 

K o 1 1 / o • 

Note here that [3, 4]-soliton has P-type interaction with other two solitons, 
which gives the last row to be (0, 0, 1, 1, 0, 0) and 3rd and 4th columns to 
be e 3 S M. 3 (see Corollary EJ- 
a5. n = (6,4,5): This case is similar to the previous one, i.e. two P-types 
and one T-type. A3 is 

'10 1\ ____. 

10-1-20, o - • 

v 1 1 1 0/ o - • 

a6. n = (6,5,4): All the vertices are of P-type, and the matrix A3 can be 
written by 

















1\ 





— 


— 


- 


- 


• 





1 








-1 


, 







— 


— 


• 










1 


1 





0/ 









• 







The Young diagram Y is maximum with three boxes. 

In those examples, one shoud note that the pattern of the zeros in the last 
three columns (non-pivot columns) presents the shape of the corresponding Young 
diagram Y~ . 

Case b: n + = (1,2,4): There are ra n + = 1-2-2 = 4 different types of 3- 
soliton solutions. Now the Y + has one box, there is one O-type interaction in those 
solutions. 

bl. n = (3, 5, 6): There are two T-type vertices and one O-type. An example 
of A3 matrix is 








-1 





1 


M 





— 


• 











1 


2 





-1 


- 2 < 







— 


— 


• 













1 


1 


1/ 











— 


• 



Notice here the duality in £(1, 2, 3) = and £(4, 5, 6) = 0. 
b2. n = (3, 6, 5): The Y~ has one box, and the vertices are all different types, 
0-, N- and T- types. The matrix A3 can be 

'10-100 2 \ o - • 

110 0-1, o - - - • 

v 1 1 / o • 

b3. n = (5, 3, 6): This is similar to the previous case of b2. A3 can be 

'l 1 2\ o - - - • 

110 0, o • 

,000111/ o - • 
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Figure 4. 3-soliton solutions. All the cases assume 



for 



all j = 1,...,6. The top left figure shows the solution in the 
case al with the parameters (fci, . . . , ke) = (—3, —2, 0, 1, 2, 3) and 
t = 7. The top right one shows the case a3 with (—3, —1, 0, 1, 2, 3) 
and t = 7. The one in the bottom left shows the case a5 with 
(—4, —1, 0, 1, 2, 3) and t = 4, and the one in the bottom right shows 
the case cl with (—2, —1,0, 1, 2, 3) and t = 4. 



b4. n = (6,3,5): The Y has two boxes, and there are two P-type vertices 
and one O-type. The matrix A3 is given by 



'1 r 

011000 

,000110, 



o — — — — 
o • 

o • 



Case c: n + = (1, 2, 5): There are m n + =1-2-1 = 2 different 3-soliton solutions. 
Since the Y + has two horizontal boxes, there are two O-type vertices in those 
solutions. 
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cl. n = (3,4,6): The Y~ has no box, and the 3-soliton solution has one 
T-type interaction with two O-types. An matrix A3 can be 








-1 


-2 


0\ 





— 


• 











1 


1 


1 


, 







— 


• 


















1 1/ 













• 



c2. n = (4,3,6): The Y has one box, and the solution has one P-type 
interaction with two O-types. The matrix A3 is given by 



1 


-1 


0\ 





- 


— 


• 






1 1 





, 







• 














1 1 













• 



Case d: n + = (1, 3, 4): The Y + has two boxes, and there are m n + = 1-1-2 = 2 
different 3-soliton solutions which have two O-type interactions. 

dl. n = (2,5,6): Since Y~ — 0, there is one T-type intersection with two 
O-types. An example of the matrix A3 is 



1 


1 











\ 





• 












1 





-1 


2 • 









— 











1 


1 


1 












d2. n = (2,6,5): Now Y has one box, and the solution has one P-type 
interaction with two O-types. The matrix A3 is given by 





1 











°\ 





• 
















1 








-1 , 









— 


— 


• 











1 


1 














• 





Case e: n + = (1, 3, 5): This is the ordinary 3-soliton solution with the matrix, 





1 

















1 


1 














1 


1 



Some of the 3-soliton solutions are illustrated in Figure 0] 

5. A^-SOLITON SOLUTIONS OF THE KdV EQUATION 

Here we consider AT-soliton solutions of the KdV equation, and show that they 
cannot have a resonant interaction. ./V-soliton solutions of the KdV equation can 
be obtained by the constraint du/dy = in the KP equation. Since line solitons 
are given by [rtf, nj]-solitons for j = 1, . . . , N, the constraint implies that all the 
solitons are parallel to the y-sxis, i.e. 

Cj = k n + + k - = for all j = 1, . . . , N. 
Now from the ordering k\ < • • • < k2N, we assume fc^'s to satisfy 

k\ < &2 < ' ' ' < ^n < 0, fcjv+j = — fcjv-j+i for j — 1, . . . , N. 
Then we take the set (n + , n~ ) as 

nf = j, nJ = 2N-j + l for j = 1, .. .,7V, 
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Figure 5. 3-soliton solutions of the KdV equation labeled by 
n + = (1,2,3) and n = (6,5,4). The left figure shows the 3- 
soliton solution in x-y plane at t = 25 with 9j = 0, Vj. The right 
one shows the solution in the x-t plane with different Oj's. The 
fc/s are given by (h, . . . , ks) = (-4/3, -1, -1/2, 1/2, 1, 4/3). 



which leads to Cj = for all j = 1, .. . , N. In terms of the matrix An, this 
corresponds to 



(1 





.4 



N 












\0 ••• 1 * ••• 0/ 

This matrix indicates that any pair of line solitons are of P-type, that is, all the 
interactions are nonresonant. Each line soliton has the following form, 



u(x, t) 



2k] sech 2 0j , 



with 9i 






k p + 6» 



Thus each soliton has the velocity dx/dt 
of the KdV equation in Figure 



-k]. 



We illustrate a 3-soliton solution 
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